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Abstract 

We investigate the stability against small deformations of strings dangling into Ac/S'5-Schwarzschild 
from a moving heavy quark- ant i- quark pair. We speculate that emission of massive string 
states may be an important part of the evolution of certain unstable configurations. 
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1 Introduction 

In El E] , string configurations in AdS^-Schwarzschild were exhibited which were argued 
to describe the propagation of a heavy-quark meson through a thermal quark-gluon-plasma. 
Briefly, the conclusion of these works is that quarkonium systems can propagate without 
experiencing any drag force, provided they are small enough and their velocity is not too 
high. Part of the interest of this topic stems from measurements at RHIC of J/ip suppression 
relative to binary scaling expectations jB] , which is not as severe as expected based on models 
assuming color screening (see for example jHj)- Could the no-drag property predicted using 
AdS/CFT have to do with the weaker than expected J/ip suppression? 
The present paper aims to address more modest questions: 

1. Are the string configurations studied in 012111] stable against linear perturbations? 

2. When / if one of these configurations is unstable, what does it evolve into? 

Along the way we will encounter an analytic expression for the shape of the string describing 
a moving quarkonium system. 

A related work [7] appeared recently which emphasizes global comparisons of different 
branches of the configuration space of strings attached to the same flavor brane. Our analysis, 
along the lines of point ^ is complementary in that stability is examined solely from the 
point of view of small deformations about a given classical solution. 

The rest of the paper is organized as follows. In section |21 we give an account of linear 
perturbations. In section El we speculate about the evolution of unstable perturbations. We 
end in section @] with a brief summary of our conclusions. 
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2 Perturbations of a quarkonium string configuration 



The metric of AdSs-Schwarzschild is 



ds 2 



z 2 



fl _ ^\ dt 2 + (dx 1 ) 2 + (dx 2 ) 2 + (da; 3 ) 2 + —^rdz 2 



We will assume that while the ends of the string are separated in the x 2 direction by coor- 
dinate distance £, they are constrained to move with velocity v in the x 1 direction. Thus we 
are already restricting ourselves to a special case of the analysis of E] , which is roughly 
the configuration of [2] without spin. It was found in [3J H] that no solutions are possible 
unless the separation between the two ends is less than some critical value I < £ max - If the 
separation is smaller than £ max , then two solutions are possible. We will argue that only 
one of them is stable against small perturbations, namely the one that dangles less far into 
AdSs-Schwarzschild and has lower energy. 

We work in the string's equilibrium rest frame: that is, we boost x 1 — > 7(x 1 + vt) and 
t — > 7(t + x l v) to get 
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The string's classical dynamics is described by the Nambu-Goto action: 
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where X^ represents the embedding of the string into Ac/S'5-Schwarzschild and a a = (r, a) 
are the worldsheet coordinates. The position of the string can then be described by 



X°(t,<t) X l (T,a) X 2 (T,a) X 3 (T,a) Z(r,a 



(4) 



We choose "static" gauge, X°(t, a) = r. An obvious completion of this gauge choice is 
X 2 (t, a) = a, but it turns out that Z(t, a) — a is a better one because then we are able to 
express the equilibrium configuration of the string in closed form (see (JJJ)). An odd feature of 
this gauge choice is that it covers only half the string: X 2 (t, a) is a double- valued function. 
In exploring perturbations we must therefore patch together two half-solutions at the string 
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midpoint. 

With the gauge choice explained in the previous paragraph, the Nambu-Goto lagrangian 
density is 
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Under the assumption of an everywhere smooth string configuration, one may demonstrate 
that X 1 = X 3 = jlj. If we solve the Z-equation of motion we get, for the left half of the 
string (where X 2 < 0), 

X» = (r X 2 Aa) a) , (6) 



with 
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where iq is the Appell hypergeometric function. The right half of the string is described by 
the same solution (J7J) with an overall sign flip, X 2 — > —X 2 . Here, X 2 (a m ) = 0, so a m can be 
understood as the largest ^-coordinate achieved by the string. As noted in (3J |3] , the above 
solution only exists as long as the string doesn't get too close to the horizon, namely as long 

as 

OVn < Z v = Z H {1 - V 2 )^ . (8) 



We have thus parameterized our solution by a m , with < a m < z v . We can relate a m to the 
coordinate distance between the string's endpoints by writing I = — 2X^(0) which gives 
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An equivalent result to (jHJ) has already appeared in jH] . A plot of I as a function of a m can be 
seen in figure d where v = 0.9. For this particular value of v the maximum £ max is attained 
at o m = (jjnax ~ Q.51zh- As we shall see, all configurations that have a rn < cr max are stable 
with respect to small perturbations, and all configurations with a m > <r max are unstable. 
To prove this claim, we do a linearized stability analysis of the equilibrium configurations 
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Figure 1: A) The separation I as a function of a m for v = 0.9. (See (jUJ).) As explained below 
(J7J), <j m is the maximum value of z along the string. The solid part of the plot corresponds 
to stable configurations; the dashed part corresponds to unstable configurations; and the 
red dot corresponds to a = <7 max , where the separation between the quark and anti-quark is 
maximized. B) A cartoon of string configurations corresponding to three points on the curve 
from A. For any separation of quark and anti-quark less than the maximum (corresponding 
to the red curve), there is a stable string configuration (solid) and an unstable one (dashed). 



given by ©. We write 



X»=(t 5X\r : a) X 2 (a) + 5X 2 (r,a) X 3 (r,a) a + 5Z(r,a)j 



(10) 



where we have set 5X 2 = because any non-vanishing 5X 2 corresponds to reparameter- 
izations of the in-plane perturbations. With this choice, and denoting 5X 1 = 5X\\ and 
5X 3 = 5X±, we derive the linearized equations of motion that follow from the action Q: 
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First note that equation (|12|) for 5Xj_ decouples and simplifies to 

-#+(l-^) fl-4)^--fl-^V)^l^ = ( 22 ) 

Callan and Giiijosa have looked at a similar equation in [Oj, where they only considered the 
pure-AdS case. It is easy to check that our equation reduces to equation (7) of [HJ, 
provided we take the pure-AdS limit zh — ► oo. 

As we shall see, it follows from equation that all transverse normal modes 5X±(t, a) = 
Re {SX^(a)e~ luJT } have cj 2 > 0, and so the equilibrium configuration is stable with respect 
to small perturbations in the x 3 direction. We can see this directly by solving the above 
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equation, but we choose a different approach: we look at how the hamiltonian of the system 
changes when we make time-independent perturbations of the shape of the string in the 
x 3 direction. More explicitly, if the hamiltonian always increases when we statically perturb 
the shape of the string, then it must be that the equilibrium configuration of the string is 
locally stable. Physically, the situation is analogous to having a particle at the minimum of 
a potential well. If we slightly change its position, the particle would tend to go back to the 
lowest potential energy configuration. 
In our case, the hamiltonian density is 
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c 7 d T x i + ^d T z-c, 



d T X l d T Z 
and, if we assume static configurations of the form 
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which is indeed positive for any small 5X±(a). Hence all normal modes in the x 3 direction 
have uj 2 > 0. The question remains whether or not the in-plane normal modes given by 
equation (fTTj) share the same feature. 

A hamiltonian approach similar to (|23|) could also be used to explore in-plane normal 
modes, but we proceed instead to solve the coupled equations in (JTTJ) directly. This allows 
us to determine the lowest normal mode explicitly. Because we're mostly interested in the 
unstable regions, we look for solutions to (jlljl of the form 



6X\\(t,(t) = Re{5X||(a)e nr } 
5Z(r,a) = Re{6Z(a)e nr } . 



(26) 
(27) 



Plugging this ansatz into we obtain the eigenvalue problem: 
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The eigenfrequencies u = iVL can be found by solving these equations and imposing appro- 
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priate boundary conditions at a = and a = o m . 

At o = 0, equations ()28|) give the following asymptotic relations for the solutions 5X\\(cr) 
and SZ(a): 

6X\\ (a) = ai - ^aifi V + a 2 a 3 + O (a 4 ) (29) 
SZ(a) = U x tf + (3 2 a + (a 2 ) , (30) 

where the choice of the four complex constants a±, a 2 , /3±, and j3 2 determines the whole series 
expansion. The boundary conditions that we want to impose are 



&Xji(0) = 



SZ(0) = 



(31) 



which imply setting a x = Pi — 0. 

The first consequence of imposing these boundary conditions is that all in-plane normal 
modes have real Q 2 . Indeed, if without loss of generality we assume that j3 2 is real, the 
reality of Q 2 follows from examining the leading order behavior of the imaginary part of 
the 5Z equation in QJ: Im{LHS} = -CV 8 Im {Qa 2 } + 0(a 9 ), whereas Im{RHS} = 
C 2 /3 2 a 3 Im {fi 2 } + 0((T 4 ), for some non-zero real constants C\ and C 2 . Since the leading 
order behavior should be the same on both sides, we conclude that Im{f2 2 } = 0, so O 2 is 
real. Moreover, if Q itself is real, as it is the case for the unstable modes, we further have that 
SXh and SZ are real everywhere, because they satisfy the differential equation which 
in this case has real coefficients. Being interested in the unstable modes, we will henceforth 
assume that all these quantities are real. 

Imposing boundary conditions at o = o m is slightly more complicated, because we have 
to patch smoothly the two halves of the string. In order to do this, we first express 6X\\ and 
SZ not in terms of our usual variable a, but rather in terms of the equilibrium coordinate 
X 2 , whose expression in terms of a is given in J7J) for the left half of the string. The patching 
conditions then read 
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(32) 
(33) 



Since we expect the eigenfunction corresponding to the lowest normal mode to be an even 



7 



function of X 2 , the boundary conditions that we impose are, for the left half of the string, 
dSXn 



dXl 



dSZ 



. (34) 



The remaining challenge is to translate the boundary conditions (jHlj) i n t° corresponding 
boundary conditions for 5X\\ and SZ expressed in terms of a. 
At a = a m , the asymptotics for 6X\\ and SZ are: 

5X\\(a) = a + aiy/a m - a + a 2 {<J m - o) + O (\a m - a\ 3/2 ) (35) 
8Z{a) = b + b iy /a m -a + b 2 (a m -a)+0(\a m - a\ 3/2 ) . (36) 

Here,the constants ao, ai, bo, and b\ can be chosen independently The chain rule then gives: 

d8X\\ da d8X\\ d5Z da d5Z , . 

(37) 



dXl dXl da dX 2 e dXl da ' 

From (J7J) we get that = C y/a m — a + (\a m — cr| 3 / 2 ) for some constant C, and plugging 
this into (pTTJ) we get 

)c ai + O {Va~^) (38) 



dXl 2 
dSZ 1 



dXl 2 



Cb x + O (v 7 ^^) , (39) 



e 



We can see from (|3*H)l and (|3*9^1 that the boundary conditions translate into requiring 
ai = &i = 0. 

In addition to the boundary conditions (j31)l and (j34j) . we should also impose a nor- 
malization condition on <JXn or 8Z, because any solution of the above eigenvalue problem 
can be multiplied by an arbitrary constant and still be a solution. We choose to impose 
5Z'(0) = 1. For numerical purposes it is useful to solve ()28|) by imposing 5X\\(0) = SZ(0) = 0, 
8X'J'(p) = Qa 2 , and 5Z'(0) = 1, and then vary the two parameters a 2 and Q until we can 
satisfy a\ — b\ — 0. 

We find that for all the configurations to the left of the maximum in figure ^ we can 
never have a± = b\ — 0. These, then, are stable configurations, for which the lowest normal 
mode has u 2 > 0. The configurations to the right of the maximum in figure however, do 
have a normal mode with uj 2 < 0, or equivalently, Q 2 > 0. 

For example, if v — 0.9 the maximum £ m ax is attained at a m ~ 0.51zjj. If we look 
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Figure 2: Eigenfunctions for v = 0.9 and <r m = 0.6zh- The derivatives 5X'^(a m ) and 5Z'(a m ) 
are finite, which means that the perturbed string configuration is smooth at a = cr m : see the 
discussion following (J34j) . 




0.52 0.54 0.56 0.58 0.62 111 

Figure 3: fl 2 as a function of a m for v = 0.9. 



at a configuration with a m = 0.6zh, we find that we can solve the eigenvalue problem 
given by (|2*Sj) with the boundary conditions (j3"Tj) and if we choose Q ~ 3.15/ zh- The 
corresponding eigenfunctions can be seen in figure 121 Moreover, the dependence of Q 2 as a 
function of a m can be seen in figure H3 The dependence of Q 2 on cr m is evidently very close 
to linear, but there are small deviations from linearity. We do not know if these deviations 
are just artifacts of imprecise numerics, or if Q 2 should actually be a linear function of a m . 

3 The evolution of unstable perturbations 

For the string configurations with a m > cr max , which we have demonstrated to be unstable 
equilibria, there are several possibilities: 

1. The string may evolve into the stable configuration with a m < cr max . 
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Figure 4: A string (purple) and anti-string (green) trailing from a heavy quark and anti- 
quark, pictured here as c and c. The string and anti-string attract, and we propose that 
configurations where they come together and self-annihilate are an important aspect of the 
dynamics. The annihilation process of the string-anti-string pair produces highly excited 
closed strings, also pictured, which then fall into the horizon. Note that the closed strings 
(light blue) may be of a different type from the string-anti-string pair: see the discussion in 
point B below. 

2. The meson may dissociate, so that its constituent quarks fly away in slightly different 
directions, each trailing a string as in [H ll()j. 

3. The string may dump energy into the horizon through a mechanism involving highly 
excited strings. 

The third possibility is the focus of this section. Briefly, we propose that configurations 
where the string and anti-string trailing from quark and anti-quark come together into an 
unstable pair may play a significant role in the dynamics. See figure |U The string-anti-string 
pair tends to self-annihilate, but depending on parameters, this may be a slow process. If 
it is sufficiently slow, then based on the numerical studies of pQ, one may expect that the 
string-anti-string pair will "blow back" into the trailing string configuration studied in [T|ll()|. 

In figureEJwe have indicated string-anti-string annihilation by letting the trailing doubled 
string break into long fragments below some elevation. This is not a process that we have 
good analytical control over. But, whether it is fast or slow compared to the relaxation into 
the trailing string shape, there are reasons to believe that the decay products are highly 
excited string states: 
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(A) If the coupling is weak, g s <C 1, then there is some small probability (of order g 2 s ) per 
unit length of the string to break. The result will naturally be long fragments of string, 
which are indeed highly excited string states. 1 

(B) If the coupling is strong, g s ^> 1, then we should pass to an S-dual picture where the 
string-anti-string pair is replaced by a Dl-anti-Dl pair. Studies of similar systems in the 
context of bosonic string theory ^2] show that in the limit where the S-dual coupling 
is small (meaning g s ^> 1 in our original language), the unstable brane configuration 
decays into highly excited open strings. Once the brane is gone, the open strings must 
close. Note that in the original language, the decay products are closed Dl-branes. 

In either the weak or strong coupling scenario, the highly excited strings (or D-strings) can 
decay into lighter states as they fall into the horizon. 

Evidently, we are dealing with considerable uncertainties on the string theory side, 
especially if we attempt to extend the picture to intermediate values of g s (recall that 
9s — Qym/^ i s essentially a s , so g s ~ 0.3, at least naively, for a comparison to the QGP at 
RHIC) . Nevertheless, it seems very likely to us that emission of massive string states plays a 
role in an AdS / CFT description of energy dissipation from mesons. Another aspect of this 
role is that if a quark and anti-quark come close enough together, their trailing strings will 
attract and then annihilate into a stable meson configuration, again resulting in the emission 
of massive string states. 

The speculative ideas discussed here could be the basis for a quantitative calculation of 

(T mn ): instead of sourcing the five- dimensional graviton with the trailing string, one should 

source it with pressureless dust, which is a reasonable approximation to the dynamics of 

massive string states. In the case of two quarks coming together to form a meson, the initial 

positions and momenta of the dust particles could be specified based on the sections of the 

trailing string-anti-string pair they originate from. 

1 In the limit of small g s , one could hope to find a classical solution where the string and anti-string 
are coincident below some elevation in AdS^-Schwarzschild and assume the standard trailing string shape 
of E3] ■ We are encouraged to learn ^I] that piecewise solutions have been found where the string and 
anti-string join at a kink. A string-anti-string pair cannot be simply joined onto this configuration because 
force balance at the junction can't be maintained unless the kind turns into a cusp. It seems likely that the 
string- anti-string attraction, weak though it is in the g s — > limit, has to play a role in the description of a 
configuration of the type we have depicted in figure 0] 
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4 Conclusions 



We have demonstrated by an explicit analysis of linear perturbations that one branch of 
Lorentzian solutions describing meson propagation through a thermal medium is locally 
stable, and that the other branch is unstable. The stable configurations are the ones where 
the string dangles less far into anti-de Sitter space: see figure ^ This stability calculation 
tends to support the overall picture presented in [21 El 15' m which it was argued that heavy 
quarkonium systems, as described in AdS/CFT, can in certain circumstances propagate 
without a drag force (at least at tree level in string theory) through the thermal medium. 

We have also speculated about the possible evolution of unstable configurations of a 
string between a heavy quark and anti-quark. The main substance of our remarks is that 
an obvious and perhaps dominant form of energy loss for such systems is the emission of 
massive string states. Such states might also play a role in an AdS/CFT description of the 
formation of mesons out of heavy quarks. 

It is essential to bear in mind that the AdS/CFT description of mesons may may or may 
not capture the dominant aspects of the dynamics of J ftp or T propagating through a real- 
world QGP. The absence of light dynamical quarks in the AdS/CFT description is the most 
worrisome contrast with real- world QCD. But even if we retreat to the most conservative 
position that AdS/CFT provides merely an analogous system to the QGP produced at RHIC, 
a more complete description of the dynamics of quark-anti-quark pairs remains an interesting 
avenue for future research. 
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